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Abstract — In most MIMO systems, the family of waterfall 
error curves, calculated at different spectral efficiencies, are 
asymptotically parallel at high SNR. In other words, most MIMO 
systems exhibit a single diversity value for all fixed rates. The 
MIMO MMSE receiver does not follow this pattern and exhibits 
a varying diversity in its family of error curves. This work 
analyzes this interesting behavior of the MMSE MIMO receiver 
and produces the MMSE MIMO diversity at all rates. The 
diversity of the quasi-static flat-fading MIMO channel consisting 
of any arbitrary number of transmit and receive antennas is 
fully characterized, showing that full spatial diversity is possible 
if and only if the rate is within a certain bound which is a 
function of the number of antennas. For other rates, the available 
diversity is fully characterized. At sufficiently low rates, the 
MMSE receiver has a diversity similar to the maximum likelihood 
receiver (maximal diversity), while at high rates it performs 
similarly to the zero-forcing receiver (minimal diversity). Linear 
receivers are also studied in the context of the MIMO multiple 
access channel (MAC). Then, the quasi-static frequency selective 
MIMO channel is analyzed under zero-padding (ZP) and cyclic- 
prefix (CP) block transmissions and MMSE reception, and lower 
and upper bounds on diversity are derived. For the special case 
of SIMO under CP, it is shown that the above-mentioned bounds 
are tight. 

Index Terms — MIMO, linear receiver, MMSE, diversity 



I. Introduction 

Linear receivers are widely used for their low complexity 
compared to maximum likelihood (ML) receivers. In the con- 
text of MIMO systems, linear receivers such as the minimum 
mean square error (MMSE) receiver are adopted in some of the 
emerging standards, e.g. IEEE 802.1 In and 802. 16e. Therefore 
the analysis of MMSE receivers is strongly motivated by both 
theoretical and practical considerations. 

A significant amount of research has focused on linear 
receivers, however, their performance is not fully understood 
in the MIMO channel. For instance, the distribution of the 
output signal-to-interference-plus-noise ratio (SINR) of the 
linear MIMO receiver is still unknown except in asymptotic 
regimes (large number of antennas, and high/low SNR) [1]- 
[4]. The outage and diversity of MMSE receiver have also 
been a subject of interest. It has been observed [5]-[7] that 
while the MMSE receiver can extract the full spatial diversity 
of the MIMO quasi-static channel at low rates, it does not 
enjoy this feature at high rates. 

Figure Q] shows the outage probabilities (for various spec- 
tral efficiencies R bps/Hz) of MMSE and ML receivers 
respectively. Clearly, one of the main differences between the 
two characteristics is the slope of the error curves, i.e., the 
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diversity. Figure Q] shows that in a 2 x 2 MIMO system the 
ML receiver achieves diversity 4 at all rates. However, the 
MMSE receiver diversity varies with the operating spectral 
efficiency. From a system design perspective, obtaining the 
MMSE diversity is important in order to understand the broad 
tradeoffs involved in the determination of the operating point 
of the system and predicting its performance. 

In this work we seek answers for the following questions: 
when can the MMSE receiver exploit the full diversity in 
MIMO channel? More generally, how does the diversity of 
the MMSE receiver vary with the system parameters such as 
spectral efficiency R, the number of antennas, and in case of 
inter-symbol interference channel (ISI), the channel memory? 

The well-known and powerful framework of diversity- 
multiplexing tradeoff (DMT) is not sufficient to answer the 
above questions, because the DMT framework cannot distin- 
guish between different spectral efficiencies that correspond 
to the same multiplexing gain. In the MIMO MMSE receiver, 
rates that correspond to the same multiplexing gain can 
produce different diversities. 

We approach the problem of MMSE reception in MIMO 
flat fading channels through a rate-dependent approximation 
of the outage probability and then proceed with bounding 
the pairwise error probability (PEP) from both sides using 
the outage. This leads to a closed-form expression for the 
diversity-rate tradeoff which reveals the relationship between 
diversity, spectral efficiency, and number of transmit and 
receive antennas. The approximation of outage and PEP as 
functions of rate requires more delicate handling compared 
with the DMT analysis, as certain ratios and terms that simply 
vanish in the DMT analysis are in our case relevant and must 
be carefully handled. 

We then analyze the frequency -selective, quasi-static MIMO 
channel. Specifically we consider single carrier (SC) MMSE 
equalization under zero-padding (ZP) and cyclic-prefix (CP) 
transmission. SC-MMSE provides an attractive alternative to 
orthogonal frequency division multiplexing (OFDM) due to its 
low complexity and natural avoidance of the peak-to-average 
power ratio problem. The use of cyclic prefix and zero padding 
has been investigated in the literature, but the explicit tradeoff 
between the spectral efficiency and diversity of MIMO SC- 
MMSE under these two schemes has been unknown and is the 
subject of our work. We show that the diversity is a function of 
number of antennas, channel memory and spectral efficiency, 
and obtain the explicit tradeoff in the special case of SIMO 
under CP transmission. 

The results of this paper fully characterize the MIMO 
MMSE diversity in the fixed rate flat quasi-static regime. 
We analyze both the cases N > M and N < M, showing 
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Fig. 1 . Outage probability of ML receiver (left) and MMSE (rigr 

that in either case it is possible for the system to be limited 
to a diversity strictly less than MN. More specifically, the 
central result of the paper is as follows: with M transmit 
and N receive antennas (for any N and M) the diversity is 
d = \ (M2~ v - (M - N) + ) + ] 2 + 1 N - M \ \ (M 2" £ - (M - 
N) + ) ], where (■)+ = max(0, •) and [•] denotes rounding 
up to the next higher integer. Our results confirm and refine 
the earlier approximate results on the diversity of MMSE 
MIMO receivers that were obtained for very high and very 
low rates [5]-[7]. The MIMO MAC channel is also studied. 

Some of the related literature is as follows. The performance 
of MMSE receiver in terms of reliability goes back to [8] 
where outage analysis was performed for MMSE SIMO di- 
versity combiner in a Rayleigh fading channel with multiple 
interferers. In the context of point-to-point MIMO systems, 
Gore et al. [9] compared the performance of MMSE D- 
BLAST with the ordered successive cancellation V-BLAST. 
They show that the former has better throughput at low- and 
moderate SNR. Onggosanusi et al. [5] studied MMSE and 
zero-forcing (ZF) MIMO receivers and noticed their distinct 
outage performance at high-SNR, specifically for large number 
of transmit antennas and low spectral efficiencies R, but 
provided no analysis. 

Hedayat and Nosratinia [6] considered the outage proba- 
bility as a function of fixed rates R under joint and separate 
spatial encoding, but for MMSE they obtained results only in 
the extremes of very high and very low rates. Kumar et al. [7] 
provided a DMT analysis for the system of [6] and observed 
that the DMT analysis does not predict the diversity of MMSE 
receivers at lower rates. We note that all existing analyses are 
limited to the case where the number of receive antennas (N) 
is greater than or equal the number of the transmit antennas 
(M). 

This paper is organized as follows. Section [XT] describes 
the system model. Section [Fll] finds the exponential order of 
outage. Section [IV] bounds the codeword error probabilities 
using the outage values, and derives the final result. Section |\H 
extends the result to the MAC channel. Section [VTl calculates 
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with M = N = 2 antennas and for rates R= 1, 4, and 10 bps/Hz 

the diversity of MIMO MMSE reception in frequency-selective 
block-transmission systems. Section [VTll provides simulations 
that illuminate our results. 

II. Linear Receivers 

The input-output system model for flat fading MIMO chan- 
nel with M transmit and N receive antennas is given by 

y = Hx + n (1) 

where H e i^nxm j s ^ c h a nnel matrix whose entries are 
independent and identically distributed complex Gaussian, 
x e C Mxl is the transmitted vector, n e C Nxl is the Gaus- 
sian noise vector. The vectors x and n are assumed indepen- 
dent. We assume a quasi-static flat fading channel and perfect 
channel state information (CSI) at the receiver (CSIR) and 
no CSI at the transmitter (CSIT), therefore transmit antennas 
operate with equal power. 

We aim to characterize the diversity gain, d(R,M,N), 
as a function of the spectral efficiency R (bits/sec/Hz) and 
the number of transmit and receive antennas. This requires a 
pairwise error probability (PEP) analysis which is not directly 
tractable. Instead, we find the exponential order of outage 
probability and then demonstrate that outage and PEP exhibit 
identical exponential orders. 

Following the notation of [10], we define the outage-type 
quantities 

P out (R, N, M) 4 P(7(x; y) < R) (2) 

d out (R, N, M) = - hm (3) 

p^oo log p 

where p is the per-stream signal-to-noise ratio (SNR). 

We say that the two functions f(p) and g(p) are exponen- 
tially equal, denoted by f(p) = g(p) when 

Km log -^) = 1^ 1 °gg(/ 3 ) 

p^oo log(p) p^oo log(p) 
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Fig. 2. MIMO system with linear MMSE receiver 



The ordering operators ^ and ^ are also defined accord- 
ingly. If f(p) = p d , we say that d is the exponential order of 



A. MMSE Equalizer 

The equalizer, denoted by W, decouples the M transmitted 
data streams at the receiver (Figure |2). The MMSE equalizer is 
obtained by minimizing the mean square error (MSE) defined 
as E[||x - W H y|| 2 ]. It is usually assumed [6], [7] that the 
number of transmit antennas M is no more than that of receive 
antennas N. In the following, we start with N > M but later 
generalize it to N < M as well. 



For N > M, using the orthogonality principle, the MMSE 
equalizer is given by [5], [11] 

W = U H (HH H + p- 1 !)- 1 

= (H^H + p- 1 /)- 1 ^ (4) 

The corresponding signal-to-interference and noise ratio 
(SINR) of the output stream k of the MMSE detector is 
1 



under either uniform or non-uniform rate assignment. 

Furthermore, it has been established [6], [7] that the zero 
forcing equalizer achieves diversity N — M + 1 under both 
joint or separate spatial encoding. 

According to Theorem [TJ a MMSE receiver operating under 
separate spatial encoding (e.g. horizontal encoding V-BLAST) 
will have no more diversity gain than ZF receiver. 

III. Outage Analysis 

We now consider the MMSE diversity where the data 
stream is first encoded then multiplexed into M sub-streams, 
each transmitted by one antenna. This approach is known to 
improve the performance compared with separate coding of 
the streams [14]. Outage occurs if the channel fails to support 
the target rate [12]. After channel equalization, the M sub- 
streams Xk are decoupled and thus the mutual information 
between the transmitted vector x and the received vector y 
given CSIR is [5] 



M 



Ik 



>kk 



1, 1 < k < M 



(5) 



>kk 



[ + pH H H 

where (-) H denotes matrix Hermitian, ( 
onal element k of the matrix inverse. 

For the case N < M, it can be shown using a technique^ 
very similar to [8, Appendix A] that the SINR expression (0 
is again valid. 

The square matrix W = H H H is random, non-negative 
definite, and obeys the Wishart Distribution [12], [13]. In this 
work, the joint distribution of the eigenvalues of this equivalent 
channel matrix opens the door to the development of our 
analysis, as is also the case in many other MIMO results. 

The equalizer output is 



7 ( x >y) = ^2 I (xk,yk) 



k=l 



denotes the diag- Thus from ® and dS), P out is given by 



Pn 



M 



^log(l + 7fc ) <R 



k=l 



Substituting MMSE SINR ( 7fc ) from © in © we get 



Pout = P[ £log(I + p>V)^ >-R 



k=l 



(8) 



(9) 



(10) 



y = WHx + Wn. 



(6) 



The signal streams of the transmit antennas may be either 
separately or jointly encoded. Separate encoding is simpler 
and has been fully analyzed [6], but we mention the central 
result for completeness. 

Theorem 1 ( [6], [7]): In a MIMO system consisting of M 
transmit and N receive antennas (N M), under separate 
spatial encoding, the MMSE receiver achieves the diversity 



d out (R,N,M)=N-M + l 



(7) 



The dependence on the diagonal elements of the random 
matrix (I + pW)^ makes further analysis intractable. We 
instead proceed to provide lower and upper bounds on the 
outage probability. In Section [IV] we will show that outage 
probability (P ou t) and pairwise error probability (PEP) exhibit 
identical exponential error. 



A. Outage Upper Bound 

Lemma 1: For an MMSE MIMO system consisting of M 
transmit and N receive antennas, under quasi-static Rayleigh 
fading, we have P out (R, M, N)<^p- d ™ t( - R ' M ' N ^ where 



1 In [8] an MMSE diversity combiner is used at the receiver in the presence 
of one transmit antenna and M interferers. 



d out {R,M,N) = 



(M2-$r - (M - N)+Y 



-i 2 



4 



\N — Ml 



(M2-™ - (M — N) + y 



(11) 



where () + denotes the max(0, •). 
Proof: 

We begin by bounding the sum in ( TTOb via Jensen's inequal- 



ity 



M M 

£ log (I + P W) ~ k l < Mlog ( £ - (I + P W)^) 



k=l 



fe=l 
1 



Mlog(— tr((l + pWT 1 )) 



M 

Mio g (— y — - 

k—1 



(12) 



where ( flZb is true because trace is equal to the sum of 
eigenvalues. 

Notice that for N < M only N eigenvalues are non-zero, 
hence ( flZb can be written as 

L f 

(13) 



where L = min(Af, N). 

Substituting ( fT3l in ( [Tol l, we have 



/ l 

Pout < P V : , , ^ M2"W - (Af - JV)- 



Define: 



pA fc 

a log A fc 
log/? 



for fc = 1, n , 



based on which we can write the exponential equality 



1 



a k <1 



(14) 



(15) 



(16) 



1 + (^1 a k > 1 

Define a = [a.\, ...,a n ] and a new random variable 

M(a) = y 1 (17) 

Q fc >l 

This definition is based on the observation that the term 



l+pA fc 



defined in ( fT6b is either zero or one at high SNR, therefore to 
characterize ^2 k 1+ L pXk at high SNR we count the ones. Thus 



y = y i+ y P a -- x 



fe=i 



Qfe>l «fc<l 



M(a) + max p 01 "^ 1 

{a k :a k <l} 



(18) 
(19) 



M(a) inherits its randomness from Ai,...,A„. The bound 
in ( fT4l i is evaluated by computing the probability of {a £ 
A}, where .4 = {a : M(a) + maxj ai . Qi<1 } p afc_1 > 
M2"« — (M — -/V) + } denotes the outage event based on 
the approximation in ( TBI . In order to evaluate the probability 
of this event we need the joint distribution of the eigenvalues, 



or equivalently the distribution of a. The distribution follows 
Wishart distribution and was initially discovered by [13] . The 
distribution of a can be easily evaluated as follows [15]. 

Let R be an m x n (to ^ n) random matrix whose entries 
are £/V(0, 1). The joint PDF of the ordered random variables 
a (defined in ( fT3T l for the eigenvalues of R H R) is given by 



»( a )=K-] n (io gP ri[p- 



(m— n+l)c*i 



n 



p _a '| 2 exp 



(20) 



where K m 1 n is a normalizing factor. 

Using the distribution of a for the defined matrix R, the 
asymptotic outage bound is 

Pout< J P(a)da 
A 



K-) n {\ogpT jf[ P - {m - n+1)ai U\P 



i.=i 



exp 



da 



p~ a J\ 2 x 



(21) 



The simplification of the integral follows from [15]. The 
term outside the integral has no effect on the exponent. The 
term \p~ ai — p~ aj \ is dominated by p~ ai at high SNR. We 
now divide the integration range into A 1 = A PI W] and its 
complement. If a. ^ A', the exponential term will dominate 
the other terms and will drive the integral to zero. If a £ A', 
the exponential term is approximately 1 at high SNR and will 
disappear. Therefore 



Pout < / Y[ P -^- n+1 ^ 

P n 

= J JJ p -(«-i+«-»)* da 



-an |2 



da 



(22) 



where 
A' 



{M (a) > M2-m - (M - N) + } 
= {en > 1, a s > 1, a s +i > 0, ...a L > 0} (23) 

and S = [(M2"w - (M - N)+) + ]. The integration region 
A' has boundaries that are parallel to nonnegative orthant K™ , 
therefore the integration over multiple variables in (l22l can be 
separated: 



Pn 



n 

i=l 



P 



-{2i-l+m-n)oLi ^ 



= p-Ef=i(2i-l+m-n) 
= /) -(S 2 + (™-n)S) 



-(S 2 + |m-n|S) 



for to ^ n 
for general m, n 



(24) 



(25) 
(26) 



= P 



which establishes the proof of Lemma Q] □ 
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B. Outage Lower Bound 

Lemma 2: For an MMSE MIMO system consisting of M 
transmit and N receive antennas (and L = min{M, N}), 
operating under quasi-static Rayleigh fading, we have 

P out {R,M,N)^p- d °^ R > M > NS > where 



d out {R, M, N) 



(M2 _ w - (M - N)+Y 



\N — Ml 



(M2—a - (M-N)+Y 



Proof: The lower bound is also based on Jensen's inequality. 
Recall 



Pout 



\ 

^log(l+ 7 fc)< J Rj 
fc=i ' 

M v 

y log r < r 

ti (i + pw)^ 1 / 
1 M 1 



R 



(27) 



Let the eigen decomposition of H H H be given by H H H = 
U H AU where U is unitary and A is a diagonal matrix that 
has the eigenvalues of the Wishart matrix W on its diagonal. 
Let the vector u^, be the column k of the matrix U and uik 
be the element I of this column, we have 

(I + P W)- fc 1 = uf(I + pA)- 1 u fe 



M i |9 



Y- 

f-r 1 



Let k = argminfc5fc. Using 
in d27) 

M 



J_y i 



f pXe 

(28) 

we can bound the sum 



1 

min fc S k 
1 



(29) 
(30) 



thus the outage bound in 
ing 



can be further bounded us- 



/ 1 M 1 
P out > P Mlog — > — r < i? 



^ P( Mlog— < R 
^k 



Sz > 2~M 



(31) 



We now bound ( T3TT > by conditioning on the event B = 
{\ u ek\ 2 ^ Xf} where a is a positive real number that is 



slightly smaller than one, i.e. a 
positive number. We then have 



1 — e, and e is a small 



S- h > 2~M > 



S k > 2--\B)¥{B) 



El 



\ u ikf 



> 2~m 



BjP(B) 
> 2-#V(B) 



1 x a 
~M f-' 1+ pA 

i M 

1 \ - a r 
— > > 2" 

M^l + pXe 

^ 1 M n\ 

y - — - > —2-- 



(32) 



M 



h y T-^r > — 2~ # - (m - jv)- 

1 + pA £ a 



(33) 



where (1321 follows because P(B) is finite and independent of 
p; this can be proved similarly to [7, Appendix A]. To make the 
upcoming expressions compact, we introduce a new variabe 

k 4 M 2 -£ - (M - AQ+ 



1 L 1 
M f-f 1 + , 



> K 



(34) 



Whenever M2 m is non-integer, the constant a can be 
chosen such that [(M2"# - (Af - A^)+) + ] = [(v 2 ~^ _ 
(A/ — A r ) + ) + ] . We note this is satisfied for all rates, with the 
exception of an isolated set of points. As long as M2~m ^ N 
we have: 



P„ 



L 1 

Y — 



pXi 



> K 



y > m 



fc=l 



pAf 



(35) 



The remaining steps follow similarly to the proof of Lemma Q] 
Thus P out ^ p- d ™t w j m ^ out j s gi ven by Lemma [2] 

On the set of isolated points M2~m g N, the right hand 
side of Eq. d35l l obeys a slightly weaker upper bound by 
replacing k with k + 1. We can combine the cases where 
M2~w is integer and non-integer to write the upper bound 
compactly as follows: 

2 

d out {R,M,N) < (M2-3 + 1 - (A/ - A r )+) + 



JV- M 



(M2-& + 1 - (M - A0 + )~ 



Inspection shows that this bound is tight against the lower 
bound everywhere except its discontinuity points. In other 
words, the upper bound is left-continuous while the lower 
bound was right-continuous at the discontinuity points. 

□ 
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IV. PEP Analysis 

Recalling that the diversity is roughly defined as the slope 
of PEP at high SNR, we now proceed to bound the PEP tightly 
from both sides using the outage results already obtained. 



A. PEP Upper Bound 

We start by a lower bound that is inspired by [15, Lemma 5] 
but requires a more careful treatment since we are analyzing 
rate, not the DMT (see the Introduction). 

Lemma 3: For a quasi-static fading MIMO channel with 
MMSE receiver we have d out (R, M, N) > d(R, M, N). 

Proof: 

Denote E for an error event, and let x <E C be the transmitted 
codeword from a codebook C of size 2 Rl where R and I are 
code rate and code length respectively. Define / = WHx that 
accounts for the combined effect of channel and equalizer. The 
transmit messages are assumed equi-probable so the entropy 
H = log \ C\ = Rl. Applying the Fano inequality [16] 

= /) > *- I tey\f = f)-Wmf = f) (36) 



Rl 



Rl 



By defining Vs for any 5 > as Vs = {/ : I(x;y|/ = 
/) < l(RS)}, and noting that H(¥(E)\f e V s ) ^ H(F{E)) 
from (|3oT l, we get 



F(E\f e v s ) > 



Rl-I(x;y\f &V S ) U(F(E)) 



Rl 



Rl 



> 6 H(P(E)) 
^ R Rl 
Also by using the definition of P out we have 



(37) 



P(/ G V S ) = P(J(x; y) < l(R - 5)) = p -*~ t (R-s,M,N) 

(38) 

For small enough values of 6 > 0, we have d out (R, M, N) = 
d ut(R—3, M, N) since d ou t(R, M, N) is left-continuous with 
respect to R. Hence, by invoking ( f3Tb and d38l , the error 
probability is given by 

F ea (R,M,N) = F{E\f e V s )F(f e V s ) + 
F(E\f i V s )F(f £ V S ) 
> F{E\f G V s )F{f e V S ) 
±(8 H(F(E)) \ 
^\R Rl P 



■ -d out 
= P 



(39) 



where we have used (i— ^^f^ 



1, which was derived 
in [10]. This establishes the proof of the PEP upper bound. □ 



B. PEP Lower Bound 

We begin by writing the error probability in terms of error 
event E and outage event O 

F m {R, M, N) = F(E\0) ■ P out + F(E, O) 



In Section IIII-AI we have shown that, based on the event 

{Efe=i T+k > M2_ * - ( M t iV ) + }' the outa s e P rob " 

ability is upper bounded by P ut^p d ° ut ■ Hence, the error 
probability can be bounded as 

F ea (R, M, N)^P(E\0) p- d °^ + F(E, O) 

<: + F(E, 6) (40) 

We intend to show that p - d °^ % F{E, O), and thus 
P er r(i?, M, N) ^ p^ d °«t which produces the following lemma. 

Lemma 4: For a quasi-static fading MIMO channel with 
MMSE receiver we have d out (R, M, N) < d(R, M, N). 

Proof: 

We begin by giving a sketch of the proof then we pro- 
ceed with the details. The first part of the proof consists 
of developing a bound on PEP conditioned on H, namely 
P[sk — > Sj\H = H\. To do this we obtain an upper bound 
of the variance of the SINR which is expressed in terms 
of the eigenvalues of the Wishart matrix W, resulting in 

P[£|H = H ] < 4cxp(-(EtiIT^) _1 )- The PEP 
is used to derive a conditional union bound on error. We 
then divide the channel events into two sets based on the 
exponential order of the eigenvalues: the set where M(et) = 
and otherwise. We apply Bayes theorem on the union bound 
using these two sets. The calculation of the terms of the 
Bayesian gives F(E, 0)^p~ AIN ^ p~ d o«* as desired. 

We now proceed in detail. We want to compute the prob- 
ability that the transmitted symbol x(k) — s; is erroneously 
detected as x(k) = Sj, 

Recalling the equalizer output given by ©, define the noise- 
plus-interference signal 



n = y - Vpx = Vp(WH - I)x + Wn 



(41) 



Using the eigen-decomposition of H and noting that 
E(n) = and E(nn H ) = I, we have 

Ma = E(h) = Vp(WH - I) = -pi (W + p- 1 !)- 1 * (42) 

R R = E(hh H ) = (W + p^ir 1 (43) 

Thus the variance of the noise sample n(k) is given by 

a?(fc) = i? fi (fc,fc)-| Mfl (fc)| 2 



(W + p-'l^-p-'iW + p- 1 !)^ (44) 



where \ph(k)\ 2 is the k th diagonal of the matrix 
E(h)E(h H ) and k counts from 1 to M. 

By defining eji = ^ J _^'| , the probability of erroneous 
detection for channel realization is given by 



F\si 



H] 



si\ 2 < \e*i(y(k) - y/ps t f 



h = h 



si\ 2 < \n k \ 2 



H = H 



(45) 

where the inequality holds since \e* l (y(k) ~ y/psi)\ ^ 

\e*My(k) - y/P»i)\ = \(y(k) - y/psi)\ = \m\- 
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Denoting the real and imaginary parts of n(k) by n r (k) ~ 
Af(p r (k),a^(k)) an d ~ N(p,i(k),af(k)) respectively, 

we then have 



(46) 

Applying the property of the Gaussian tail function Q(x) Sj 
e (-x /2) f or tne pairwise error probability, we obtain 



(^| aj - 3 ,|-^r-( fc » 



s£ e 



(^| aj - ai |+Mr-(fc)) 2 



(^g|s 3 --s,|+ Mi (fc)) 2 



s£ e 



,r 2 7k} 



C-#l»i-«jl-MiW) 



i(«>) 



(47) 



where the last step holds as cr^(fc) = <7j?(fe) + cr?(fc) ^ 

Now we proceed by showing that pj(fc)^ps. Consider the 
eigen decomposition of 

[W + p -1 !] -1 = U H [A + p- 1 !]-^ 

= U g [diag{ 1 }]U (48) 

where U is unitary matrix, and A is the eigen decomposition 
of W. Note that A fc + or Afc + -i <p. Therefore, 

all elements of the matrix ±U ff [A + p _1 I] _1 U, being linear 
combination of { \ k + p -i }■ cannot grow faster than 0{p), and 
thus the elements of ±ps [W+/0 _1 I] _1 cannot grow faster than 
O(pi), i.e. ±p„(fc)^ pi and therefore p5±p„(fc) = pi. The 
same result holds for p r (k) and pi{k). 

As a result, for any Sj and Si, ^\sj — s;| ± p r (fc) = pi ± 
/Lt r (fc) = p 2 " and similarly *£-\sj — sj| ± Pi(fc) = pi. Thus 
from (|47| |. we have 



'[a* -> Sj|H = if] < 4e~ 



(49) 



Now we bound the variance in (1441 and apply it in d49l 



< E r 

L 



(W + p-'l^-p-'iW + p- 1 !)^ 



E 

fe=i 



1 + p\ k (1 + pA fe ) 2 



P 2 A fc 



Denoting the error event E and using 
of erroneous detection in (|49l is bounded as 



P[£|H = ff] sC 4e 



. 5Zfc = l (14 



g(l + pA fc)2 
(50) 

the probability 



(51) 



Applying the union bound, we get 



¥{E\H = H)<^2 



Rl \ £fc=l 



(52) 



Based on ( 1521 ). we can evaluate P(E, O) in d40b as follows. 
Recalling the exponential inequality 



/C— 1 Ct fc >l CKfc<l 



M(a) + max p° 

{Q fe :Qfc<l} 

M(a) 



(53) 



(54) 



Consider the two regions: {a : M(a) — 0} and {a. : 
M(a) ^ 1}. At high SNR the event O is equivalent to {a : 

M(ol)<^\M2-tt -(M-N)+]}. 

In the first region {M(a) = 0}, at any rate R we 
have {a : - (M - N) + ] %M(ot) = so there is no 

outage. 

In the second region {M(a) 1} the exponent order of the 
outage probability depends on the rate. We investigate these 
two regions separately. 

In the region {a : M(a) = 0}, we have max^ a,}. < 1 
since all a' k s < 1, From d52l and (|54| i we conclude that 



"(£;,0|M(a) = 0)^2 m e- py 
= 2 Rl e'^ 



(55) 



Since exponential function dominates all polynomials and 

1 — maxfc ctk > 0, we get 



lim ■ 

p— too 

which in turn yields 



-MN 



= 



\E,0\M{a) = 0)<:2 m e- p 



(l-„, a x fc c h ) 



-MN 



(56) 



We next show that the same result holds for the other region 
{a : M(a) ^ 1}. 

Following the same line of argument as we did for (l56l l but 

for M(oc) 1, we have 



P(E,0\M(a) > 1) < 2 m e 



C 9 Rl p V^-fc=i (i+ P A fc ) 



e 



s Efc 1 + pA 



^1 since M(ck)^1 

t 1 



< e e i*r(<*)p" 



(57) 



x 



^p 



-MN 



(58) 



where dSTT i is direct application of ((54) for M(a) ^ 1, 
and (|58j follows from the fact that |1 - a k \ > 1. Note that 
d58l is true for any code length I. Invoking the results of I 
and d58l ). we can now evaluate F(E, C) as follows 



F(E,0) = J ¥{E,0\M{a) = 0)P(a)da 

M(a)=0 

P(E,0\M(a) ^ l)P(a)da 

Af(a)>l 

^p- MAr J P( a )d a + p- MN 
M(a)=0 



M(a)>l 



-MN 



(59) 

J (a)da 

(60) 
(61) 



Therefore, P(J5, 0)-^p MN for all regions of a. Finally, 
( |40b becomes 

P erT (i?,Af,iV) < P(£|0) /9" dout + P(£,(5) 



= P out (R,M,N) 



(62) 



which establishes the lemma. □ 

From Lemma [3] and Lemma [4] we thus get 

Theorem 2: For MMSE MIMO Receiver under quasi-static 
channel and joint spatial encoding, the pairwise error proba- 
bility (PEP) and the outage probability P out are exponentially 
equal and the diversity gain is d(R, M, N) = d OM t(R, M, N), 
where d out (R, M, N) is given in dTTT >- 

V. Multiple-Access Channel (MAC) 

We now extend the result to the MAC channel. Consider a 
MIMO MAC channel with K users, M transmit antennas per 
user, N receive antennas (there is no condition on M, N and 
k). Assume flat fading MIMO channel, the system model is 
given by 



K 
i=l 



n 



H e X 



(63) 



where H; € C NxM is the user i channel matrix whose 
entries are independent and identically distributed complex 
Gaussian, H e = [H1H2 . . . Hrt] is the overall equivalent 
channel matrix, Xj £ C Mxl is the transmitted vector of user 



i, X = [x{ Xj . . -xjf] is the overall transmitted vector, and 
n G C Nxl is the Gaussian noise vector. The vectors X and 
n are assumed independent. We keep the same assumptions 
about the channel. That is we assume a quasi-static flat fading 
channel and perfect CSIR and no CSIT We have the following 
theorem 

Theorem 3: In a MIMO MAC system with MMSE receiver 
consisting of K users, M transmit antennas per user and N 



receive antennas, the lower and upper bounds on the per user 
diversity are respectively given by d™ AC (R) and d^ AC (R), 



d r c (R) 



[M2- R / M - (M- N) + Y 



+ 



N - KM 



{M2- R ' M - (Af - N)+Y 



(64) 



d MAC 



(R) = 



(KM2~ R / KM - (M-N)+)~ 



N -KM 



{KM2- R ' KM - (M - N)+y 



(65) 

From d64l i it is straightforward to verify the single user 
case. The machinery of the proof is mostly similar to the 
single user case. However, the outage upper and lower bounds 
are obtained in a different manner that is pointed out in the 
following analysis for N M. The case N < M can be 
similarly obtained. 

A. MAC Outage Upper Bound 

The user i outage probability can be written as 

iM 



J2 Iog(l + 7 jt)<i2 . (66) 



fc=(i-l)M+l 

where ji is the SINR of the stream k of user i. Specializing 
this to MMSE receiver we get 

iM 



l°g(I + pH e H H e )^ >-R). (67) 



k=(»— 1)M+1 

Using Jensen's Inequality the outage probability can be 
bounded as 



iM 

Pl ut ^ P(l0g (J] -(I + pHe^He)^) > 



fc=(i-l)M+l 
KM 



M ' 



^ P ( 1 °g(E^( I + ^ H e HH e)^)>^|) (68) 



k=l 



KM 

= P(Y >M2~v) 



(69) 



where ( I68t is true since the summation in the left-hand side 
of the inequality adds more positive terms (recall that (I + 
pH e H H e ) is a positive definite matrix [12]). Following similar 
steps that were used to obtain (l26*T l we can easily show that 
Pout ^ P~ dh - where d™ AC is given by d64l i. 



B. MAC Outage Lower Bound 

The outage probability can be lower bounded as follows 



pi _ 

out 



iM 



£ lDgtl + pH^He)-, 1 >-R 

fe=(i-l)M+l J 
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KM 



> ¥ (J2 l0g(I + pHe H He)^ > -R 



k=l 
KM 



(70) 



^p(^(I + pH e «H e )^>^2^) (71) 

^ k=l ' 

where ( |70l > is a trivial bound based on dedicating all KM 
antennas to one user, and dTTl l uses the same technique as in 
Section ITH-B1 and a is a positive number slightly less than one. 
Following similar steps that were used to obtain (l26l we can 
easily show that P l . 



_ JM 



where d^\ AC is given by d65l l. 



VI. Frequency-Selective Channel 

Broadband wireless systems usually operate in frequency- 
selective channels where, in addition to the spatial diversity 
obtained in MIMO broadband systems, frequency diversity can 
be achieved. Broadband systems usually employ orthogonal 
frequency division multiplexing (OFDM) or single carrier (SC) 
transmission [17]. Specifically, SC was shown to be attractive 
for broadband wireless channels due to its lower complexity, 
lower peak-to-average power ratio and reduced sensitivity to 
carrier frequency errors compared to OFDM [17], [18]. 

In this section, we investigate the diversity achieved by SC- 
MMSE receivers for two block transmission schemes, namely 
cyclic prefix (CP) and zero-padding (ZP) schemes. The CP 
and ZP are commonly used for guard intervals in block quasi- 
static channels. Although CP was initially proposed for both 
single carrier and multi-carrier systems, ZP was lately shown 
to be an attractive alternative for both systems [19], [20]. 



A. System Model 

We consider a general MIMO system in a rich scattering 
quasi-static environment. The equivalent baseband channel is 
given by multipath model with v paths referred to as the 
ISI channel in the sequel. The (y + l)-tap channel impulse 
response between the transmit antenna m and receive antenna 
n is denoted by the vector h TO „ = [/i„ m ,o, h mnA , . . . , h mn ,u\- 
We assume a block-fading model where h m „ remains un- 
changed during a transmission block. Assuming M transmit 
and N receive antennas, the received vector y& at time instant 
k is given by [10], [21] 



Yk = 



Xfc- 



n fe 



(72) 



where is the M X N channel matrix that has h mn ,i as 
its (m, n) element, Xk-i is M x 1 transmitted vector at time 
index k — i, yk is the N x 1 received vector and n^, is the 
N x 1 Gaussian noise vector at time index k. 

Consider a transmission of Ld + L e spatial vectors each of 
size M X 1, where L d is an integer representing the number of 
transmissions over the quasi-static channel and L e is the length 
of data extension to avoid inter-block interference, in the form 
of either zero-padding or cyclic prefix. The receiver discards 
the L e vectors in the case of cyclic-prefix transmission [21]. 



Stacking the transmitted vector in an M (Ld + L e ) x 1 vector, 
we can write the stacked M(Ld+L e ) x 1 transmitted as follows 



■L s )+L d +L e -V 



x k - l X k(L d +L e ) 

We can then rewrite (|72l as 

y cp = Hx + n (73) 

where y cp is the NLd x 1 received vector, x is the M(Ld + 
L e ) x 1 transmitted vector, n is the white Gaussian noise vector 
G C NLdXl and H is the channel matrix given by 

Hq Hi ■ ■ • H,, 



H 











Hi 
Ho 



Hi 




H (/ 



Ho Hi 






H, 



(74) 



The linear data extension operation maps the data vector x 
to the transmitted vector x and is shown by 



x = U. 



where U cp is given by 
U™ = 



f-MLd 

0, 



(75) 



(76) 



l ML e ¥ML e x(L d -L e )M 

The system model in d73l can now be written in terms of the 
unpadded data vector x and an equivalent channel matrix H e 
as follows 

y cp = H e x + n (77) 

where in a CP system, H e = HU cp is a NLd x MLd block 
circulant matrix constructed by block circulations of the matrix 
[H o ,H 1 ,...,H„0,...,0] T 
For the zero-padding transmission, we can rewrite (|72l as 

H e x + n (78) 



' zp 



where y zp is the N(Ld + L e ) x 1 received vector, x is the 
MLd x 1 transmitted vector, n is the white Gaussian noise 
vector G C N ^ Ld+Le ' xl and H is the channel matrix given by 



Hp 



Ho 

H, 







Hi 
H„ 







H . (79) 



H, 



Assuming perfect channel state information at the receiver 
(CSIR) and that the channel remains unchanged during the 
transmission of Ld + L e vectors, the MMSE equalizer W is 
applied to decouple the received streams (after removing the 
L e extension vectors in case of cyclic -prefix transmission). 
The MMSE equalizer is given by 



W= (p^I + He^He)- 1 ^" 



(80) 



and the unbiased decision-point SINRs of the equalizers output 
for detecting the k th transmitted stream are 

7fc= 7Z TT-tTTZ — - 1 k = l,...,ML d . (81) 



(I + pH e w H 
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In the following sections we analyze the outage diversity 
for the ZP and CP systems. The PEP analysis follows in a 
direct manner as in the flat fading case so we omit it. 



B. The Zero Padding MMSE Receiver 

It is known that in a point-to-point single-antenna ISI chan- 
nel, linear receivers can achieve full multipath diversity under 
zero-padding transmission [20], [22], [23]. In this section we 
investigate the similar question for MIMO systems whose 
receivers use linear MMSE operations in both the spatial and 
temporal dimensions. We provide lower and upper bounds on 
diversity. The bounds are not always tight, but the diversity is 
fully characterized for SIMO systems. 

We begin by analyzing the tradeoff between the spectral 
efficiency R and the diversity of MMSE receiver in the 
single-antenna ISI channel d\ s [MSE under ZP transmission. 
Tajer et al [10] shows that d^^SE var i es with R under 
CP transmission and MMSE equalization, in particular, for 
a quasi-static single-antenna ISI channel with v + 1 taps, the 
diversity of the SC-MMSE receiver under CP transmission is 
d*MMSE = l+ mm (^ \ %~ R Ld\), where L d is the transmission 
data block length. We show that the same is not true for ZP 
transmission. 

Lemma 5: For a quasi-static single-antenna ISI channel 
with v+1 taps, the diversity of the SC-MMSE receiver under 
ZP transmission is dff MSE = v + 1 irrespective of R. 

Proof: See Appendix lAl ■ 

We proceed with lower and upper bounds on diversity for 
MIMO ISI channel. 

1 ) Diversity Upper Bound: Applying the MMSE equalizer 
given by (l80l to the received vector in (|T7K the effective 
mutual information between x and Wy is equal to the sum 
of mutual information of their components [5] 



ML d 



J(x,Wy) 



1 



ML d 

E 

k=l 



Thus the outage probability is given by 
/ i MLd 

Pout 



— £ log(l+ 7fc ) <i?) 
d fe=i ' 

, ML d 

>^ lOK Tj T 

L d £> ^(I + pHe^He)-, 1 



< R 



(82) 



(83) 



ML 



1=1 



Ulk\ 



pXe 



Let k = argmirifc S k - we can bound the sum in 

ML d -, ML d 

1 ^ — ^ 1 1 ^ — > 1 

ML d ^ (I + pHe^He)-, 1 = JIT d ^ i S l 



1 



min fc S k 
1 

thus the outage bound in (l84l l can be further bounded 



(85) 



Pmi.t ^ 



Mlog 



1 M 1 



> P Mlog — < R 



1 



< R 



S k > 2-^1 



We now bound ( 1861 ) by conditioning on the event 



M 



I = ML r] - M + 1, 



ML. 



(86) 



(87) 



where a is a positive real number that is slightly smaller 
than one a = 1 — e±, and ei is a small positive number. We 
then have 



Sk > 



S- k > 2-b\B P(B) 



ML d 



Y- 

4^ l 



> 2'm 



B)F(B) 



\Utk\ 



ML d 

=ML d -M+l T H 
ML d 



> 



B)F(B) (88) 



E 

'L d -1 
ML d 

E 



e=ML d -M+l 
MLd 



£=ML d -M+l 
MLd 



l + pX e 



El Mr 
> _2-w 
1 + pX e a 



> 2~m p(g) 



> 2~m 



=ML d -M+l 



(89) 



(90) 



Mlog > 77 T < R (84) 

where we have used Jensen's inequality as in Section llH-BI Let 
the eigen decomposition of H e H H e be given by H e ff H e = 
\J H AU where U is unitary and A is a diagonal matrix that has 

the eigenvalues of the matrix H e H H e on its diagonal. Let the Note that H e "H e is not a Wishart matrix, hence the 



where (188 1 follows by removing some of the elements of the 
sum corresponding to the largest eigenvalues. The steps used 
to obtain Eq. d89l are similar to the steps used in Section lTH-BI 



eigenvalues of H e ff He be given by {Xe} with Ai ^ A2 • • • ^ 
XML d - Let the vector be the column k of the matrix U, 
we have 

(I + pH e ff H e )^ = Uf (I + pA)- 1 ^ 



analysis of Section [TT] does not directly apply here. The block 



diagonal elements of H e H H e are similar and are given by 

V 

D = ^HfH l . (91) 
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The matrix H e ff H e is Toeplitz and Hermitian. Moreover, 
the matrix D given by (f9T~b is a Wishart matrix[l 

Observe that the probability in ( |90t depends on the M 
smallest eigenvalues. We now bound these eigenvalues with 
the eigenvalues of the matrix D via the Sturmian separation 
theorem [24, P. 1077]. 

Theorem 4: (Sturmian Separation Theorem) Let {A r ,r = 
1, 2, . . .} be a sequence of symmetric r xr matrices such that 
each A r is a submatrix of A r+ i. Then if {Afe(A r ) , k = 
1, . . . , r} denote the ordered eigenvalues of each matrix A r 
in descending order, we have 

A fe+ i(A i+ i) < Afc(Aj) < A fc (A i+1 ). 

For our purposes, we consider a special case of the 
Sturmian Theorem by constructing a set of matrices 

A M , Am+i, • ■ • , A Ld M starting by the largest one A LdM = 
H e fl H e and making all other matrices Aj to be (successively 
embedded) i x i principal submatrices of H e ff H e , such that 
the smallest matrix is Am = Dj, d . Then we repeatedly apply 
the first inequality in the Sturmian to get: 



1 ML d v 
— ^l0g(I + P H o ff H e )^>-i?) 

d fc=l ' 



M log ■ 



1 



ML d 

a k 

fe=l 

ML d 



E 

i 



(I + P H e H He)" 1 > -R 



(94) 



(I + P U e H U 



>kk 



> -R 



1 



fi 1 

M 

k=l n 



P x k (n e H H e ) 
i 



P x k (B. e H H e ) 



M 



k=l 



1 

p\ k (n e H n e 



> M2" 



L d M - M > M2" 



> M2~ 



(ML a 



(95) 

M) 



\ML d (A AfLd 
^ML d -l{^ML d 

>^ML d -M+l(A.ML d 



< AjvfL, 



,-l(A M i (i -l) 
,-2(A M L £i -l) 



< 
< 



< Am (Am 

< Am 



(96) 

where d94t follows from Jensen's inequality and ((95} follows 
from setting the smallest L^M — M eigenvalues to zero. 

Now we repeatedly use the second inequality in the Stur- 



mian theorem to get 

i(Am) 

< ■ ■ 



< AA/ Ld -M(A AfL(i _i)< • • ■ < Ax (A M ) 



Am (Am) 
Am-i(A-m) 



< 



< A 
< A 



M—l (AmL^ 



< A M (A Mi J 

< AM-l(AML d 



This implies that the smallest M eigenvalues of H e ff H e 
are bounded above by the M eigenvalues of D, respectively. 
Hence: 



Ai(A M ) < 

A 



with A.ML d = H e ff He and Am = D, similar to the earlier 
case. Therefore the largest M eigenvalues of H e H H e are 
bounded below by the M eigenvalues of D, respectively. 
Therefore 



< Ai(AML d 

A 



-l) < Ai(A M lJ 



M 



^ 1 

fc=l 



1 



- > — 2' 
pX k (D) a 



(92) 



D is a sum of {y + 1) central Wishart matrices each with N 
degrees of freedom and with identity covariance matrix, i.e. 
D e W((u + 1)N,I). Therefore the analysis of Section M 
applies here and we have the following lemma. 

Lemma 6: In a MIMO quasi-static frequency-selective sys- 
tem (with channel memory v) consisting of M transmit 
and N receive antennas, the MMSE receiver diversity under 
joint spatial encoding and zero-padding transmission is upper 
bounded as 

2 



Mlog 



M 



1 



M ^ 1 

fe=i 



1 



pA fc (D) 



> Q 



(97) 



<M2~ 



+ i-{M-N)+y 



where Q = max(0, M2~^ - (ML d -M)). Recall that D is 
a Wishart matrix, therefore the analysis of Section |n] follows 
and we obtain the following lemma. 

Lemma 7: In a MIMO quasi-static frequency-selective sys- 
tem (with channel memory v) consisting of M transmit and 
N receive antennas, the MMSE receiver diversity is lower 
bounded as 

I 2 



+ \(v + l)N-M\ 



M2- 



+ 1 — (m - N)+y 



\QY + \(u + l)N-M\\Q\ 



(98) 



(93) 

2) Diversity Lower Bound: We can upper bound the outage 
probability as follows. 



Pn 



1 



ML d 

E 

k=l 



log(l + 7fc ) <R 



2 Let W(n, 53) denote a Wishart distribution with degree of freedom n and 
covariance (also called scale) matrix A n Y °f me diagonal block matrices 
Dj given by 19 It follows a Wishart distribution since if Bi S W(ni, 53) 
and B 2 G W(n 2 , E) th ™ B i + B 2 e W(m + n 2 ,E)- ^ 



under joint spatial encoding and zero-padding transmission. 

Q = max (0, M2-tt - (ML d - M)). 

Remark 1: Notice that both lower and upper bounds differ 
only in the second term of Q, i.e. (MLd — M). The diversity 
lower bound for L d = 1 is tight against the upper bound, but 
for L d > 1 the lower bound d98l > is trivial. 



C. The Cyclic Prefix MMSE Receiver 

For the single-antenna ISI channel under CP transmission, 
the explicit tradeoff between spectral efficiency and diversity 
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DeMux 
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Exten. 



H 















MMSE 




MUX 









SISO Decoder 



x(L d +L a ) x(L d +1) x(L d ) x(1) 
Extension Data 

Fig. 3. Single-carrier block transmission in a frequency-selective channel. In the case of CP, the extension is removed at the receiver prior to equalization. 



was found [10] to be d^f MSE = 1 + min(V, Y2~ R Ld\). In 
this section, we extend the analysis to the MIMO case. The 
system model is shown in Figure [3] We start with the general 
M x N MIMO system. 

The system model is again given by ( TTTb where H e = 
HU cp and x is generated by taking the IDFT of the informa- 
tion vector x [25], i.e. 



X Qtx* 



(99) 



where Q Tx is the augmented DFT matrix given by Q Tx — 
Q®Im, where 1m is the identity matrix, Q is the normalized 
DFT matrix, and (g> is the Kroenecker product. 

The NLd x MLd block-circulant matrix H e has eigen 
decomposition H e = Q Rx AQ Tx , where Q Rx = Q ® In- 
Both Q Tx and Q Rx are unitary matrices. The block diagonal 
matrix A is given by 



/Bi 



A = 



\ 



B 2 



(100) 



V B L J 
where the matrix Bj. is given by [26] 



for k = 1, . . . , Ld 



(101) 



i=0 



and Hi is the instantaneous MIMO channel (cf. Section [VI-A| ). 

Analogous to the proof of [10], we first consider the case 
where the transmission data-block length is equal to the 
number of channel taps, i.e. L4 = v + 1. In this case, the 
entries of Bfe's are i.i.d. normal complex Gaussian. 

1) Outage upper bound: The outage probability of the 
MMSE receiver is given by 

, ML d . 

,^g' 0g( (. + >>H."H.a ><R ' ll02 > 



P — 
1 mil 



k=l 
ML d 



■ (I + P H e 
-^log((I + P H e H H e )- fc 1 )>- J R 

d k=l ' 

( MLd 1 \ 

^ PfMlog ^(I + pHe^He)-, 1 > -R) 

V k=l d 



(103) 



k=l ' 

= P( ^tr^ + pBfB,)- 1 >ML d 2-Si 



i=l 
L d M 



(1 + C 



i=l k=l 



> ML d 2~M 



(105) 



Where ( |103b follows from Jensen's inequality, ( |104t follows 
from the eigen decomposition of H e , and Afc^ is fc-th eigen- 
value of the i-th Wishart matrix B^B^. 

Recall from Section [HI] that the eigenvalues of a Wishart 
matrix have the asymptotic property 

Et^Ei +!>'-' («*> 

k=l y K a k >l a k <l 

based on which we established in Lemmas [T] and |2] the 
following 

M 

where is defined in (fTBT l and s, M, and A are arbitrary 
integers. Define 

Qfc,i>l 

are i.i.d. discrete random variables with the following 
asymptotic distribution (cf. Section [HI] Equations (l22b-([26l)) 

P(f3i = n t ) = p -(«?+|JV-M|n«) for n< = 1, . . . , M (108) 

Using d!07t , the outage probability in ( |105t can be evaluated 

as 

• L d M 

P. . 1 I VV7: — > ML d 2~H 



^^(l + P A fc ,) 
id 



(109) 



where = |~M £^2 *r ] . Evaluating the probability in (11091 1 
in a combinatorial manner, we get 



fe=i 
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p-«+\N-M\nx) _ _ p -{nl+\N-M\ nLd ) (110 ) 
-(n?+|JV-M|ni) _ _ _ p -« d + \N-M\n Ld ) 



ni,n2,...,nL 



: max 

ni,n 2 ,...,n L 



(111) 



where S [0, M] for (z = 1,2,..., L<j) is the value of the 
i-th discrete random variable 0j, and (II 1 11 1 is true since the 
summation in ( II 10b is dominated by the maximum element. 

Let the set {n* k , k = 1, . . . , £d} be the set of indices of the 
optimal solution of (II 1 11 1. The set {n* k } is obtained by solving 
the following optimization problem 

- ^(n 2 k + \N-M\n k ) 



ni ,n 2 ,...,n L 



fe=l 



subject to nfc = SI 

fc=i 



< n fc < M 



or equivalently, 



min n? 

d k=l 



(112) 



subject to rife = SI 



n fc > 



The problem in dl 12l i is a quadratic integer-programming 
(QIP) problem (see e.g. [27] ). Integer programming problems 
are in general NP-hard. However, due to the simple structure 
of the objective function in (II 121 >. we can efficiently solve it, 
thus obtain a closed form expression for {n^} and hence (111 lb . 

Lemma 8: For the QIP given by (111 2b . the optimum solu- 
tion is given by: 

n* = u for 1 i ^ t 

n* = u + 1 for t + 1 < j < L d 

where u — and t — L d {u + 1) — SI. 

Proof: See Appendix iBl ■ 

Using Lemma [8] we can now evaluate the outage upper 
bound given by (111 11 1 as 



(113) 



where d cp = Q(2u + 1) - uL d {u + 1) + \N - M|fi and 

2) Outage lower bound: The bound is obtained using the 
same steps to obtain the lower bound in Section IVI-B. II It can 
be shown that 

/ 1 MLd \ 
Pout = Y, ^Sd 1 + PHe^He)^ 1 ) > -Rj (114) 



L d M 



> ML d 2~M 



(115) 



The bound in dl 15b is the same as the upper bound in dl05l l, 
thus the bound is tight and the diversity is given by dl 13b . The 
PEP analysis follows in a manner similar to Section |IV] 

Recall that so far we have considered data block length 
L d = v + 1. It can be shown that the diversity for any L d > 
v + 1 is upper bounded by the computed diversity for the case 
L d = v + 1. This bounding is derived from (1104b via FFT 
arguments similar to those used in [10], which we omit for 
brevity. A tight diversity lower bound for data block lengths 
L d > v + 1 remains an open problem, except for the SIMO 
system as discussed in the next section. 

3 ) Diversity of CP Transmission in the SIMO Channel: 

Theorem 5: In a SIMO quasi-static frequency-selective 
channel with memory v, N receive antennas and data-block 
length L d , the MMSE receiver diversity is d^MSE = 
N vci\\v(v + 1, [2~ R L d \ + 1) under joint spatial encoding and 
cyclic prefix transmission. 

In order to prove Theorem |5] we first analyze the case of 
L d = v + 1 and then generalize the result for L d > v + 1. The 
system model is given by ( |77] i where the NL d x L d equivalent 
channel matrix is given by 



h 





hi 
h„ 



hi 







hi h 2 



(116) 



where hi (for i = 0, 1, . . . , v) is N x 1 SIMO channel. Note 
that the diagonal elements of (H e ff H e ) are identical and 
equal to Ya=q h f h, . Thus the MMSE SINR for each output 
information stream is 



1 



7fe = 



(I + H e H H e ) 



- 1 = 



1 



>kk 



- 1 
(117) 



Evaluating the outage probability as in ( 11021 



i ^ d 

a k=l 
1 L " 

log— y 

L d ^ 

k=l 



1 



(I' 



P n e H n 



< R 



e >kk 



(I + P H e H H 



1 



i (l + pA fc ) 



>L d 2~ 



< R 



'kk 



(118) 



(119) 



where (II 18b follows from (II 17b and ( 11 19b follows similarly 
to JT05Y 



In a manner similar to (1105b we have = B^Bfc because 
now B is simply a N x 1 vector. For the case L d = v + 
1, the eigenvalues {A^} are distributed according to Gamma 
distribution with shape parameter N and scale parameter 1, 
i.e. Afe ~ r(JV, 1). For L d > v + 1 the Gaussian variables 
in Bfe are no longer independent and thus analyzing this case 
requires the unknown distribution {Afc}. Instead, we indirectly 
show that the diversity of L d = v+1 also holds for L d > v+1. 

Lemma 9: In a SIMO quasi-static frequency-selective chan- 
nel with memory v, N receive antennas and data-block length 
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L d = v + 1, the MMSE receiver diversity is d ( lf MSE = 
N([L d 2~ R \ +1) under joint spatial encoding and cyclic 
prefix transmission. 

Proof: 

The outage probability can be written as 



HM{a) > L d 2~ R ) 



-R\ 



(120) 



where we use M(a) = J2a k >i 1 fr° m (1 1 061 > . We thus need 
to evaluate P(a > 1). The probability density function of 
is 



fx h (x) 



1 



■x N ~ 1 e- x 



T{N) 

The distribution of a k is thus given by 

The cumulative distribution function of at is 

^a fc (x) = / /«* (J/) dy 



(121) 



(122) 



1 



r(iv) J p 
l 



iW Wo 

AT-l 



JV— 1 -rj 

r e ar 



N—l -r j 

e ar 



- :c k 



(123) 



iV-l -r j 

e ar 



it! 



(124) 
(125) 



where we have made a change of variables r = p~ x in (1123b . 
and evaluate the integral according to [24, P.334 and P.336]. 
Thus we have 



N-l 



P(a k >l) = l-e-"Y J 



P 

fc! 

fc=0 

1 —N\ 



1-1 rP 



(126) 
(127) 



where (1126t follows from the Taylor expansion for (1125) . 

From the independence of {A^}, and subsequently the 
independence of {oik}, we conclude that M(a) in d 1 20t is 
binomially distributed with parameter p~ N . Hence, similar 
to [10], we have 



a i 



>L d 2~ 



\M(a) > L d T 



F(M(a 

4 =LL d 2-«J+l 
L d 



p-"*(l-p-») 



N\n—i 



^ p -N({L d 2- R \+l)_ 
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Fig. 4. Outage probability of MMSE Receiver, M = N = 3 for R=l, 1.5, 
2, 3, 4.5, 4.8, 5, 10 bps/Hz 



which concludes the proof for L d = v + 1 ■ 

For L d > v + 1 we follow steps similar to [10]. 

Lemma 10: Consider two SIMO systems both operating 
under quasi-static frequency-selective channels with memory 
v. One system has data block length L dl > v + 1 and the 
other L d . 2 > L dl , we have the following property 



fc=i v 



1 



(l + pA fc ) 



> m \ — 



fc=i v 



1 



(1 + pAfc) 



> TO 



for any m£K. 

Proof: The proof has similarities with the SISO case 
developed in [10, Lemma 2], but is not a trivial extension (see 
Appendix let. ■ 

Using Lemma [10] and the results in [10, Theorem 2], 
Theorem [5] is established. 

VII. Simulation Results 

Simulations generate Monte Carlo random channel realiza- 
tions and calculate outage probability by checking the appro- 
priate linear MIMO receiver mutual information for the quasi- 
static flat fading model. Figure |4] shows the case M = N = 3. 
According to Theorem [2] d out — 1 for R 4.755, d out = 4 
for 4.755 > R ^ 1.755, and d out = 9 for R < 1.7549. 
Figure |4] shows the diversity step between R — 4.5 and 
4.8bps/Hz. The slope of diversity 9 is difficult to measure 
precisely with simulations, but it is approximately observed. 
Figure shows the outage probability for R — 1,4 and 10 
with the lensen bound, with a diversity transition at R = 2. 
Figure [6] shows the case of M = 2, and N = 3 again with 
transition at R = 2. In Figure [7] simulations results for N = 2 
and M = 3 are given and compared with N = 3 and M = 2. 
Theorem [2] gives the diversity for both systems. It is observed 
that when N > M the break point of the slopes occurs before 
its counterparts in M > N case. Lower rates were difficult to 
simulate precisely. 
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Fig. 5. Outage probability of MMSE Receiver, M = N = 2 for R (left to 
right)= 1,4, 10 bps/Hz 
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Fig. 6. Outage probability of MMSE Receiver, M = 2, and N = 3 for R 
(left to right)= 1.5, 2.5, 4 bps/Hz 



R(left to right)=1 .8, 4, 10 bps/Hz 
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Fig. 7. Outage probability of MMSE Receiver for both cases N > M (solid) 
and M > N (dashed). The spectral efficiency R (left to right)= 1.8, 4, and 
10 bps/Hz 



VIII. Conclusion 

This paper settles the long standing problem of the diversity 
of the MMSE MIMO receivers under all fixed rates for any 
number of transmit (M) and receive (N) antennas, giving the 
result as d = \M2~v - (M - N)+] 2 + \N - M\ \M2~v - 
(M — N) + ~}. The analysis confirms the earlier approximate 
results [6], [7] showing that the system diversity can be as high 
as M N for low spectral efficiency and as low as N—M +1 for 
high spectral efficiency. The result is extended to the multiple 
access channel (MAC). We also analyze the case of frequency- 
selective MIMO channel under cyclic-prefix and zero-padding 
transmission, and obtain the explicit tradeoff between rate and 
diversity. 



Appendix 



A. Proof of Lemma \5\ 



Consider a single-antenna ISI channel h = [ho,...,h v ], 
where v is channel memory. The transmitter sends a block of 
Ld + v symbols (i.e. the extension L e = v), the last v symbols 
of which are zeros to remove the inter-block interference. The 
system model is given by 



y = H e x 



n 



(128) 



where x is the transmitted length- (Ld + v) vector. We consider 
the case where the padding length is equal to the memory of 
the channel. The results are also valid for L e > v as a direct 
result of [10, Theorem 2]. 

The outage probability of MMSE receiver under ZP trans- 
mission is given by [10] 



Pn 



1 



— y io g (— 



i 



pH e H H 



T~i)< R ) 



e/kk 



^i> (1+ (H.»H.)- 



P 



< R) 



k=l 
1 



e Jkk 



! > -R) 



L r ,2 



-R 



^tr(H e ff H e 



< 9 



(129) 



(130) 



(131) 



where (|1291 l represents the outage probability of zero-forcing 
equalizer which upper bounds that of the MMSE. The bound 
in ( 1130b follows from Jensen's inequality. 

We want to show that tr(H e ff H e ) _1 in ( 1131b is propor- 
tional to ||h||~ 2 . Thus it is straightforward to obtain full- 
diversity at any R since [15] 



< 9 



--P 



-La 



(132) 



where c is a constant that is independent of h. 

To show that this is indeed the case, we use the result 
of Tepedelenlioglu [22], [28] which provides a family of 
linear zero-forcing equalizers that is capable of achieving 
full multipath diversity in zero-padded systems under certain 
constraints. We paraphrase the result for convenience. 
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Lemma 11 ( [22], [28]): Under zero-padded transmission, 
there exists a family of left-inverses of H e , denoted by G, 
such that ||G||~ 2 ^ C||h|| 2 for some constant C independent 
of the channel vector h. Moreover, we have ||W^jr|| ^ ||G||, 
for any G satisfying GH e = I, and W^f is given by 

W ZF = (He^HJ-'He^. (133) 

Applying the ZF equalizer Wzf on the channel output 
given by ( 1128b we get the equalized signal y = x + z, where 
z = Wzjfii. The filtered noise power P z can be evaluated as 

P z =Etr[zz H ] 

= tr[E((H e ff H e )- 1 H e H nn H H e (H e H H e )- 1 )] 

= trffHe^He)- 1 ] (134) 

where we assume the noise is uncorrelated and has variance 
equal to one. 

Using the properties of the Frobenius norm, P z can be 
bounded as 



P,=E(||W»/n| 



< E(||W 



|n|| 2 ) = L d ||W z/ | 



(135) 



Using ( 1134b . d 1 35b and Lemma [TT] the trace in (1131b can 
be bounded by 



trKHe^He)- 1 ] ^L d ||W 



zf\ 



C\\h\\ 2 ' 
Thus from ( 1131b we have 

Pom s=;p(c 2 ||h|| 2 < P ~ l ) 



(136) 



= p-^ +1 \ 



(137) 



where Ci = C 2 R is a constant independent of h and p. 

Note that the constraints and construction methods in [22], 
[28] for the zero-forcing equalizers to achieve full multipath 
diversity in ZP systems do not apply in CP systems. That is, 
Lemma [TT] is not true for CP transmission. This is because 
the equivalent channel in CP systems does not have the same 
properties that were used in [22], [28]. 

B. Proof of Lemma \Ej[( QIP Problem) 

Consider the following Quadratic Integer Programming 
(QIP) problem 



mm 

ni ,ri2 , ■ --^l 



E»2 

fc=l 



(138) 



subject to n k = 

fe=i 

n k > 0. 



9. 



where 51 and I are integers. 

Consider a candidate solution vector [rii, . . . , nk, ■ ■ • , ng\. 
We partition the variables in this vector according to their 
values into + 1 sets JVj = {n^ : n k — j} for < 
j < ft; clearly some of these sets may be empty. Denote 
the membership of each set Sj = \Afj\. Furthermore, let 



SI = ml + K where m is the divisor and K is the remainder 
of the division of O by I. From the constraint in (I138t we 
have 

t u n 

n k = J^jSj =m£ + - m)Sj =m£ + K. (139) 

k=l 3=0 3=0 

Evaluating the objective function: 

t o 

J2 n k = E( m + j~ m ) 2 ^ 

fc=l 3=0 

a n 
= £m 2 + 2m V^(j — m)Sj + — rn) 2 Sj 

3=0 3=0 
Q 

= lm 2 + 2mK + ^2{j-m) 2 Sj (140) 

3=0 
Q 

^ Im 2 + 2mK + -m)Sj (141) 

3=0 



= lm 2 + 2mK + K 



(142) 



where (TT401 and (TT421 use E"=o(j ~ m ) S 3 = K < which 
follows from (1 1 39b. 



We now propose that one may achieve optimality when all 
variables take values either to or to + 1. In that case, 

n k = mS m + (m + 1)(£ - S m ) =m£+(£ — S m ) 

k 

n l = m2 S m + (m + l) 2 (£ - S m ) = em 2 + 2mK + K. 

k 

where we substituted the value of l—S m from the first equation 
into the second equation above. This shows that the variables 
taking values morm+1 achieves the lower bound in ( 1142b . 
At optimality S m = (m + 1)1 — Vl. 

C. Proof of Lemma [77J 

We begin by showing that for any integer multiplier of 
Ld x = v + 1 denoted by Ld 2 = TLd 1 (T £ N) and any 
real- valued to £ (0,2^ ), we have 



L dl 

E — - 

q= 



^(l + pA g ) 



> TO \ — 



1 



> TO 



(143) 

Note that for SIMO-CP system, \ q = h q h q , where b q is 
the TV x 1 vector given by 

b« =Y,h n e~ 3 ^T forg=l,...,L di (144) 

n=0 

where h„ is the channel gain as a function of the tap index n, 
and the superscript i = 1, 2 is used to distinguish the variables 
in two systems with data block lengths L dl and Ld 2 . 

Recall that we can take a -point signal and apply a L d2 - 
point DFT on it (after zero-padding), which will result in a 
resampling in the Fourier domain at Ld 2 points. Following [10] 
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we can write the explicit relationship between entries of t/ 1 ) 
and t>( 2 ) as 



»=i 

where 



V =E & ifM 9 = 1,2. ..,£ d2 and/ = l,2,...,iV. 

(145) 

. (g-l)2,rl, dl 



1 1 - e ^ 

dl 1 - e ^ ^2 ^1 ' 



a. 



(2) 



(i) = '"gl^'.ll 2 
q,l logp 



Define a)!', = '"V" 9 '" ■ Note that bil) ^ , = &[, 2) and 



Q o 1 i > for 3 = M- since = ■ 



"T(g-l),J 



T(q-lU-T "8.! 

From (1 145b , we have 

iCi 2 - E w a i*2 ) i a +EE^K? • < 146 > 

i—l i— 1 s— 1 



We now analyze each part of the sum in (11461 ). For the set 
of indices A = {i : i = T(k - 1) + 1, Jfe = 1, . . . ,L d *.}, the 
coefficients {V^} are non-zero constants, then l^il 2 !^ ; 1 2 = 
V/. Noting that r\ must be real- valued, and defining 
a v = Eq. < fT46b can be written as 



',(2) 1 2 



t=l 1 " 

M 



(147) 



Note that if r\ < the second term in ( 11471 i should be 
smaller than the first term since otherwise the right-hand side 
of ( 11471 ) will be negative while the left-hand side is positive. 
Thus for r\ < we have a„ mirij a,^- Also, for a^Owe 
have p~ min * q m + A min - °£j . Thus we always have 

p-***"™ + iLp-av^p-vto*"™ , leading to the following 
lemma. 

Lemma 12: For ct^j and offl defined above we have: 
p~ a fl %p~ mia * °# < min, ag 3 for q e A. 

We now partition the DFT points into two sets A — {T(i — 
1) + 1,* = 1, . . . ,L dl } and B = {1, . . . , £d 2 }\{T(i - 1) + 
1 , z = 1 , . . . , } We now define the event: 



A 



V = {mina^ < 1 , mina^ < 1 , ... , mina^ < 1} 

i ' i ' i ' 

and proceed to evaluate the probability 



„(!) 



where (11481 ) follows since X q = h^h q and Si and S% are 
given by 



* = E 



i 



— ' 1 i 1-a , 

^ = E at rz^w 



(1) 



We now evaluate (11491 

S1 + S2 >m) =P( Si + S 2 >m 



V x P(D) 



Si + S2 > m 



Note that subject to the event T>, we have 

1 



V ) x P(X>) 

(150) 



* 2 = E 



iV 1-c 



qeB l + L/=iP 
Therefore this term can be asymptotically ignored. Also sub- 
ject to T>, we have 



ft=E- 



1 



9=il + E;=iP 9 ' ; 
and since with probability one, L dl > m, the other (non- 
negative) term can be asymptotically ignored. Thus, both the 
terms involving the set B can be altogether ignored and we 
have: 

y n~~r^ > TO ) =^(si> m\v)F(D)+ 



Mi 



Si > m\V)F(V) 



§(1 + P\) 



> m 



We have thus established ( I143l l when L^L^, We must 
now show that the same result holds for any T' when L dl \ T'. 
To do so, let Ld 2 = T'L dl , then we have 



1 



1 (1 + PA 9 ) 



> m = 



y 1 



> m 



(151) 

Using ( 1143b when L dl and (11511 ) when T'|L ( j 2 together 
establishes ( 1143b for any two positive integers. 



^(l + A,) 



> m 



E 



ti + pE^I&Sl 2 



-,N 1,(1)19 + 2^ 



^i + pEz=iI^I 2 ^i + pELI© 2 



= P 5i + 5 2 >m 



> m 
(148) 

> m 

(149) 
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